Evidence for spin memory in the electron phase coherence in graphene 
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We measure the dependence of the conductivity of graphene as a function of magnetic field, tem- 
perature and carrier density and discover a saturation of the dephasing length at low temperatures 
that we ascribe to spin memory effects. Values of the spin coherence length up to eight microns are 
found to scale with the mean free path. We consider different origins of this effect and suggest that 
it is controlled by resonant states that act as magnetic-like defects. By varying the level of disorder, 
we demonstrate that the spin coherence length can be tuned over an order of magnitude. 
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In a diffusive conductor at low temperature, the weak 
localization (WL) correction L] to the conductivity orig- 
inates from the scattering of a coherent electron around 
a closed loop, Fig.[lja). The electron wave traverses this 
loop in both clockwise and ant i- clockwise directions and 
interferes with itself at the point of intercept. If con- 
structive, the wave is localized at the intercept and has 
less probability to contribute to the current. A weak 
magnetic field will affect the phase accumulated by the 
electron wave, which tends to destroy the interference 
effect and result in positive magnetoconductivity (MC). 

In graphene, the WL effect has proven to be a powerful 
tool in understanding electron scattering and dephasing 
0-11] ■ Graphene has been predicted to be an ideal spin- 
tronic material because of its low intrinsic spin-orbit and 
hyperfine interactions and shown to be the first 

material to achieve gate-tunable spin transport 11- 15j|. 
However, previously measured spin lifetimes are orders 
of magnitude shorter than expected from the intrinsic 
spin-orbit interaction 12, 1^ 17 1. In a manner similar 
to inelastic dephasing, loss of spin memory can limit the 
size of the trajectories that contribute to the WL effect. 
In this Letter, we show how spin memory effects can be 
resolved in the WL effect through measurements of the 
conductivity as a function of magnetic field and temper- 
ature at different carrier densities. 

Our conductivity measurements in graphene devices. 
Fig. [ijb), were performed in a constant current regime 
in the temperature range from 0.02 to 5K. The car- 
rier density, n, was controlled by applying a voltage, 
Vg, to the n-Silicon gate: n — ^Vg, where the con- 
stant 7 = 7.7 X 10^" cm"2/v was determined from the 
Shubnikov-de Haas oscillations. Figure [2{a) shows the 
conductivity, a, as a function of at a temperature of 
30 mK. An almost linear dependence of ct on n was ob- 
served over the entire range of n for which the WL correc- 
tion was measured. This correction was determined from 
measurements of the low-field MC and the temperature 
dependence of the conductivity. 

In the measurements of MC, the magnetic field was 
varied over a 20 mT range at a fixed carrier density. (To 
suppress the mesoscopic conductance fluctuations that 
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FIG. 1: (Color online) (a) A coherent electron wave shown 
scattering off impurities to form a closed trajectory, (b) 
Schematic of our graphene samples, where the Si02 thick- 
ness is 280 nm. The table shows parameters of three studied 
samples, where dimensions are given in fim and mobilities in 
cm^/Vs. The mobilities are shown for three regions of the 
carrier density [see Fig. [2l[a)]. 



result from the flnite sample size Jl8|. the resistance at 
each magnetic field was averaged jS!, l4| over several close 
values of n in each region.) The MC of sample SI at 
different temperatures is shown in Fig.[5Kb). Positive MC 
is observed at all carrier densities for the whole range of 
studied temperatures. Analysis of the MC [Sj, |4| using 
a theoretical model developed for WL in graphene y 
allows us to determine the effective dephasing length L^ps- 
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Here Sa{B) = a{B) ~ cr(0), F{z) = Inz + V(0.5 + z'^), 
ij{x) is the digamma function, r^^ = 4eDB/h, D is the 
diffusion coefficient, t^^ = + , r^p and Ts are the 
electron dephasing and spin coherence times, and 
are the elastic inter- and intra-valley scattering times. 
The observed positive MC in all samples indicates a WL 
correction to the conductivity. There is no hint in our 
results of a transition to weak antilocalization that could 
result from strong spin-orbit interaction l^, 2^ or weak 
inter- and intra-valley scattering Q . 
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FIG. 2: (Color online) (a) The dependence cr(n) measured at 
30 mK in sample SI. Red circles represent regions where the 
MC was measured. The inset shows a single Dirac cone occu- 
pied up to the Fermi level, £f > ksT. (b) The MC in SI in 
region III. (c-e) The temperature dependence of the dephas- 
ing length in SI extracted from the MC. Different symbols 
represent different experimental runs. The solid lines are fits 
to Eq. [2] at high temperature. The dotted lines show L^^* . 



The results of the analysis using Eq. [T] are shown in 
Fig. WLs)-{e). One can see that above ^ 1 K the dephas- 
ing length, L^s = \JDt^s, decreases as the temperature 
increases and agrees with the theory (solid line) of de- 
phasing due to electron-electron interaction in the diffu- 
sive regime 
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where ep is the Fermi energy, r is the momentum relax- 
ation time and a is a prefactor. This regime corresponds 
to the condition ksTr/h < 1, which means that two 
interacting electrons experience many collisions with im- 
purities during the time of their interaction H/ksT. This 
inequality is always fulfilled under the studied experi- 
mental conditions. We used Eq. [2]to fit the temperature 
dependence of L^s {T) at T > 1 K and found that a ss 1 
at all carrier densities, as expected from theory [2l| . 

At temperatures below 1 K, the values of the dephasing 
length deviate from theory (Eq. [5]) and eventually satu- 
rate. The value of the saturated dephasing length, L^^^ , 
increases with carrier density. This saturation results 
from a temperature-independent contribution to L^s(r), 
which we attribute to electron spin effects, characterized 
by the spin coherence length Lg. To draw this conclusion, 
we must also consider three other potential causes of such 
saturation: (i) sample size (ii) electron overheating 



[18l |. and (iii) coupling of microwave radiation [22|. (i) 
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FIG. 3: The temperature dependence of the conductivity in 
samples SI (a-c) and S2 (d-e) at different carrier densities. 
The solid curves are fits to &a = &a^^ -f &a^^^ . Different 
symbols in (d-e) correspond to different experimental runs. 



contacts of the sample, and, therefore, the maximum pos- 
sible dephasing length is the sample length, L. Since 
U^l ^ L this mechanism can be simply dismissed, (ii) 
Electrons can be overheated by a high source-drain cur- 
rent. In our experiments we have verified that the 1 nA 
current used at the lowest temperatures does not cause 
overheating since we see a temperature dependence of 
the conductivity over the full experimental range down 
to 20 mK [see Fig. We will consider the effect of (iii) 
later when we analyze the dependence of L^^l on n. 

The WL correction also manifests itself in the tem- 
perature dependence of the conductivity, (5(t(T), at zero 
magnetic field, which is shown in Fig. |31 The mea- 
sured conductivity, (5cr(T) = a[T) - a{To) = 1/ p{T) - 
1/p{Tq), where Tq is the lowest studied temperature, 
decreases logarithmically down to 20 mK. Solid lines in 
Fig. El are fits to 5a = Jct^l _^ ^^eei^ analyze 
the WL correction, we first explicitly determined the 
Altshuler-Aronov (electron-electron) interaction correc- 
tion, 5a^^^{T), which arises from the scattering of elec- 
trons from a random pattern of Friedel oscillations. We 
did this by two independent methods, as described in the 
Supplementary Material: measu ring the temperature de- 
pendence of the Hall coefficient [23[ , and by suppressing 
the WL correction with an applied perpendicular mag- 
netic field 0). Thus, there are no fitting parameters 
related to the contribution of (5(t™^ in Fig. [3l From 
the MC analysis using Eq. [T] we determined the inter- 
valley and intra- valley scattering lengths. We found that 
Li ^ 1 /im and <^ I /im, which are both much smaller 



Dephasing of electrons can occur in the source and drain 



than L;p. As a result, they have a negligible effect on the 
form of Sa{T). To fit the contribution of calcu- 
lated values of at temperatures below 1 K using Eq.O 
where the prefactor a is determined from the higher tem- 
perature (T > 1 K) dependence of L^s{T), Fig. [2I^c)-(e). 
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FIG. 4: (Color online) The spin coherence length as a function 
of the mean free path determined from the MC (circles (SI), 
squares (S2) and hexagons (S3)) and (t(T) (diamonds (SI) and 
triangles (S2)). The solid line is a linear fit to the data. The 
inset shows the dependence of Ls on the inter- valley scattering 
length for sample SI in regions I, II and III. The dashed line 
is a guide to the eye. 



Thus, the only fitting parameter is Lg, which appears 
from Fig. [3] to scale with carrier density. 

To investigate the effect of disorder, we compared val- 
ues of Lg with those of the elastic mean free path Ip. 
As typical values of Ip in graphene flakes tend to be 
~ 100 nm, we increased the range of Ip by weakly irradi- 
ating sample S3 with gallium ions to introduce additional 
scattering sources [25| whilst maintaining diffusive trans- 
port. As a result, our full range of Ip extended over an 
order of magnitude from 20 to 180 nm. Values of Lg ex- 
tracted from the analysis of MC and o'(r) are shown in 
Fig. 21 as a function of Ip. The values determined from 
the MC and a(T) studies in different samples and ther- 
mal cycles are consistent and all lie on the same line 
that increases linearly with Ip. These data were taken 
at different n ranging from 0.2 to 5.5x10"'^^ cm~^, so it 
is clear from the figure that there is no dependence of 
Ls/lp on n. This rules out possibility (iii) above, that 
saturation of L^ps could be caused by the coupling of mi- 



crowave radiation to the weakly localized electrons [22 1 



We estimate that this would result in L^^^/lp oc l/\/n 
at high frequency, or oc ^/n at low frequency, which is 
clearly inconsistent with our experimental results. 

We now consider the effect of spin memory in graphene. 
In semiconductors, the presence of spin-orbit interaction 
(SOI) will tend to produce an antilocalization effect [l9| . 
As the temperature is lowered, this would manifest itself 
as a change in the sign of the MC as increases be- 
yond the relevant SOI length scale jl^l- However, this is 
not observed in the current experiment. An alternative 
scenario, unique to graphene, is that SOI conserves the 







Ip due to 




-^s / 


weak 


strong 


long-range 




short-range 


short-range 


Coulomb 


intrinsic SOI 




\/n 


\/n 


short-range SOI 


1 


1/VH 





TABLE I: Approximate dependence on carrier density n of 
the ratio of the spin-orbit relaxation length to the mean free 
path Ls/lp using estimates from recent literature [tI. [27l - [3H . 



z-component of electronic spin, Sz such that 

SH = g^zS, + ^,[S„y„(r) + ApUp{r)] . 

The first term is the Kane-Mele intrinsic SOI with 
g ^ 1 /ieV , which is far too small to explain our obser- 
vations. The second term is a short-range perturbation 
formed by C-C bond distortion caused, for example, by 
adsorbents [26[. Parameters Va, a = x,y, z, describe the 
strength of coupling with sublattice 'isospin' as indicated 
by 4x4 Hermitian matrices Eq, whereas Up, P = x,y^, 
describe coupling to valley 'pseudospin' matrices Ap [2]. 
On the one hand, these terms lead to random preces- 
sion of the electron spin around the z axis, causing the 
relaxation of the in-plane polarization of electrons. On 
the other hand, they distinguish between up and down 
spin electrons causing an effect identical to time-reversal 
symmetry breaking for Dirac electrons. Their effect is 



equivalent to a substitution — > t^^ + in Eq. [l] 
where the rate describes the relaxation of the x-y 
component of spin. 

To assess the plausibility of the idea that the satura- 
tion of the observed dephasing rate is due to SOI, we 
compare, in Table U the predicted density dependence of 
the elastic mean free path Ip = vpT with that of the SO 
relaxation length Lg for various assumptions about the 
nature of disorder and type of SOI, taking into account 
that Ls/lp oc Ts/t. The almost-linear dependence of 
the measured conductivity on carrier density [Fig. HJa)] 
indicates that the momentum-relaxation rate increases as 
the density is lowered, r^^ oc l/^/n. This is consistent 



with the presence of remote Coulomb scatterers |28l - l30j . 
or strong short-ran ge s catterers such as resonant states 
or lattice vacancies 127. 311 



Moreover, it indicates that 

momentum relaxation is not dominated by weak short- 
range correlated disorder (for which r""'^ oc -^/n). This 
rules out the possibility that the observed saturation is 
due to SOI because, as shown in Table U the observed 
linear relation between Lg and Ip would require the dom- 
inance of a weak, short-range correlated disorder, both 
in elastic and SOI channels. 

In metals, saturation of the dephasing length is often 
attributed to spin-flip scattering due to the presence of 
magnetic atoms [s^ . Since magnetic atoms possess an 
internal degree of freedom (their spin) which fluctuates 



4 



with time, collisions with them will tend to cause phase 
and spin relaxation on length scale Lg. We attribute the 
observed saturation of in our graphene samples to the 
same spin-flip scattering from 'magnetic' defects. In the 
context of graphene, such defects may be chemisorbed 
magnetic atoms, hydroxyl groups, vacancies, or resonant 
localized states in the vicinity of the Dirac point 33- 
IstI ]. The consistency of our results for different samples 
and thermal cycles allows us to dismiss the influence of 
magnetic atoms. Also, in the presence of such atoms we 
would not expect the dependence Tg oc ej? to hold true for 
all samples. In contrast, the presence of resonant states 
(which can result, for example, from hydroxyl groups and 
vacancies) can explain our observations. For a resonant 
state at |eo| ^ ep and such small localization radius tq 
that e^/ro ^ sp, the amplitude of a spin- flip process 
involving the exchange of an electron between the con- 
duction band and the localized state is Afi oc l/ep- This 
would lead to the spin relaxation rate oc Sp^ which 
would result in Lg/lp independent of carrier density. 

To investigate further, we examined the behavior of 
the inter-valley scattering length Li which is expected to 
arise from atomically sharp defects. For sample SI, we 
found that Li decreases with increasing carrier density, 
as shown in the bottom inset in Fig. HI At the same 
time, Ls increases (in agreement with Tg cx ep)- For 
the ion-irradiated sample S3, we observed a pronounced 
defect-induced peak in the Raman spectrum (see Supple- 
mentary Material) and a significant decrease of Li with 
respect to the non-irradiated samples indicating that ion 
bombardment produced vacancies. Nevertheless, S3 ex- 
hibits the same correlation between Lg and Li as exists 
in samples SI and S2 demonstrating that vacancies can 
act as 'magnetic' defects in graphene. 

In conclusion, we have measured the weak localization 
correction to the conductivity in graphene and found that 
the dephasing rate at low temperatures is limited by the 
spin memory. We have demonstrated that the spin co- 
herence length can be tuned experimentally by showing 
that it has a direct relation to the mean free path. By 
comparing with predictions of the density dependence 
of is and Ip, we suggest that spin decoherence is domi- 
nated by spin-flip processes caused by resonant states at 
the Dirac point. It was shown that these states can be 
caused by vacancies in the crystal and act like magnetic 
defects. Our values of Lg of up to ~ 8 fim clearly show the 
promise of graphene for future spintronic applications. 
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